Majorana edge states in interacting one-dimensional systems 
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We show that one-dimensional electron systems in proximity of a superconductor that support 
Majorana edge states are extremely susceptible to electron-electron interactions. Strong interactions 
generically destroy the induced superconducting gap that stabilizes the Majorana edge states. For 
weak interactions, the renormalization of the gap is nonuniversal and allows for a regime, in which 
the Majorana edge states persist. We present strategies how this regime can be reached. 

PACS numbers; 



Introduction. The possibility of realizing Majorana 
bound states at the ends of one-dimensional (ID) con- 
ductors formed by topological insulator edge states, semi- 
conductor nanowires or carbon nanotubes in the prox- 
imity of a superconductor |l|-|8|, as well as by quasi-one- 
dimensional superconductors [9| has led recently to much 
activity. An important factor for the interest is the po- 
tential application of the Majorana edge states as ele- 
mentary components of a topological quantum computer 
0, ■ In a nanowire the Majorana edge modes exist 

because of the p-wave nature of the induced supercon- 
ductivity, which is the result of the projection of the su- 
perconducting order parameter onto the band structure 
of the wire, consisting of helical, i.e., spin (or Kramers 
doublet) filtered left and right moving conducting modes. 
In such a setup, the Majorana edge states appear as 
particle-hole symmetric Andrcev bound states at both 
ends of the wire, with a localization length ^ inversely 
proportional to the induced superconducting gap A, and 
their wave function overlap is typically proportional to 
exp(— L/^) with L the wire length. The independence 
and the particle-hole symmetry of the two bound states 
is only guaranteed if this overlap is vanishingly small, 
therefore large L and A are required. 

Electron-electron interactions renormalize the proper- 
ties of one-dimensional conductor and so modify A as 
well as the localization length of the bound states. In this 
paper, we study these interaction effects in systems with 
helical conduction states that are in contact with a su- 
perconductor. We show that superconductivity and Ma- 
jorana edge states are stable only at weak interactions. 
Strong and long-ranged interactions generically suppress 
superconductivity and so delocalize and suppress the Ma- 
jorana edge states. For weaker and screened interac- 
tions, superconductivity and the Majorana edge states 
remain stable only if the renormalization flow reaches 
the strong coupling limit for the induced superconduct- 
ing gap at a correlation length ^ ^ L. This regime is 
reached for a large induced superconducting gap, best 
possible screened interactions, and the longest possible 
wire length L, which outlines the necessary strategy in 
the experimental search for Majorana edge states. Un- 



der these conditions, although the electron interactions 
in most cases substantially reduce the size of the gap, the 
Majorana edge states remain strongly localized at each 
end. 

In the following, we first illustrate the effect of elec- 
tron interactions on the Majorana bound states using 
the fermion chain model of Rcf. [10]. In particular, we 
show that for strong interactions the gap can entirely 
close and the system becomes equivalent to a gapless free 
electron gas. Motivated by this insight, we turn to a con- 
tinuum theory for the nanowires, allowing us to include 
the interactions more effectively and to move beyond the 
restriction to a half-filled chain. 

Fermionic chain. The prototype model for Majorana 
edge states is a one-dimensional open lattice of sites i = 
l,...,N described by the model ^ [H] 
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where q arc tight-binding operators of spinlcss fcrmions, 
for example the electron operators of the helical conduc- 
tion bands, i > is the hopping integral, A > the 
triplet superconducting gap, ^ the chemical potential, 



and 



4. 



In terms of the Majorana fermion ba- 



sis 
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(t ± A)/2. 



I), the model is 



At t = A and 



rewritten as H 

«f E^=l7^^7^^ with W± 

/i = 0, the only nonzero interaction is and the ground 
state corresponds to pairing of Majorana fermions be- 
tween neighboring sites 7i^7i+i, with an excitation gap of 
2u)4-. In the open chain, 7^ and 7^ no longer appear in 
H and remain unpaired. They form the two Majorana 
bound states that are localized on a single lattice site at 
each edge of the wire and can be occupied at no energy 
cost. For /X 7^ or A 7^ i, the two edge Majorana modes 
are coupled to the bulk system and their spatial exten- 
sion becomes larger, on the order of ^ a/ In 
with w = max{|/i|, and a the lattice constant. In 

the finite system, the overlap of the two Majorana states 
at both ends of the chain is proportional to e^^""^^, and 
the two states are independent only for Na S> ^. 
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In such a system, interactions between the fermions 
critically affect the existence and stability of the Ma- 
jorana edge states. Indeed, they lead not only to a 
further coupling of the Majorana edge states to the 
bulk system, but also can substantially reduce the 
bulk gap size. As an illustration, we include into the 
model the repulsive nearest neighbor interaction H' = 
UY.'^Ji {ni - 1/2) (rii+i - 1/2), with U > Q. The mean 
field contribution of this interaction is inessential. The 
direct part can be removed by further tuning ^ as well 
as the edge potentials, while the exchange part adds to 
w± and the contribution can be removed by tuning 
t — A. Quantum fluctuations, however, cannot be sup- 
pressed, and their role reaches much further. 

Indeed, it is straightforward to show that interac- 
tions can entirely close the superconducting gap. For 
strongly interacting t — A = C//4 we can map H 
by a Jordan- Wigner transformation to the spin chain 
H = tY.^Jii.(^f(^t+i + where erf ^'^ are spin 

1/2 operators (normalized to ±1) defined by a = |(crf + 
'io'i)Y[j<i '^j ■ ^ further Jordan- Wigner transforma- 
tion to new fermion operators c; = ^(crf 4- icrf ) Y[j<i '^j 
we then see that H = —'2tJ2^=i^{clci+i + c-^^^Ci), which 
describes a free gapless fermion gas in which the localized 
states have disappeared. 

The role of interactions are therefore crucial for under- 
standing the stability and existence of the Majorana edge 
states. In the following we use a continuum description 
for a quantitative analysis, which allows us to include the 
interactions more effectively, first at half filling as in the 
discrete model, then away from half filling. 

Continuum model. For the continuum theory, we fo- 
cus on a quantum wire with Rashba spin-orbit interac- 
tion in a magnetic field with proximity induced singlet 
superconductivity I^-IgJ . The non- interacting part of the 
Hamiltonian for the quantum wire can be written as a 



sum of two parts, Hq = Hq^'' + i^o^^ where Hq'^' is given 
by (throughout the paper h ~ 1) 



(2) 



r(i) 



(2) 

where is the electron operator for spin a, the summa- 
tion over repeated spin indices, a,/3, is assumed, r is the 
coordinate along the wire, p = —idr, is the spin-orbit 
velocity, and Az is the Zeeman energy of the magnetic 
field applied along the spin z direction perpendicular to 
the spin-orbit selected spin x direction. The second part, 
, includes the induced singlet superconducting term 
with order parameter As and is expressed as, Hq = 
i J drAs^l,alfj'^i/2 + h.c. Without interactions, i7, 



(1) 




has the eigenvalues e± = /2m ± {urpY + {Az/2Y 
and corresponding eigenmodes 5'±(p). Expanding the 
singlet superconducting term in this eigenbasis leads to 
superconducting order parameters of the triplet (within 



and subbands) as well as of the singlet type (mix- 
ing and ^1/+ subbands). The Majorana edge states 
require triplet pairing 0-l2, 15, which is achieved by 
tuning the chemical potential to lie within the magnetic 
field gap such that only the subband is occupied. In 
Ref. [6[, Majorana edge modes were derived using the 



full Hamiltonian H, 



(1) 



("2") 

m^> and were shown to exist in 



the limit Az > \/ Ag -I- \j? . The same physics is also ob- 
tained by restricting to the occupied subband, which 
will be assumed in the following. For Az ^ As, ankp, 
with kp ~ VmAz, the pairing then takes the compact 



form US, [11, m 



« {Alkp) / dr*L(r)p«'L(r) +h.c. 



(3) 



with the effective triplet superconducting gap A 
As{oiRkFlAz). 



In the following we work in the diagonal basis jl7| with 
the fermions confined in the r > region. The open 
boundary condition forces the fermion fields to vanish at 
the boundaries, thus ^_(r = 0) = ^'_(r = L) = 0, where 
L is the length of the wire. In terms of the slowly vary- 
ing right, 7?.(r), and left, Or), moving fields, the field 
^'-(r) acquires the form, ^'_(r) = sin(fcr)c_ (fc) = 



n{r) 



'ikpr 



C{r), where c-{k) is the annihilation 



operator in the XP- subband. We note that TZ{r) = 
—C{—r). Thus, the kinetic energy can be expressed in 
terms of TZ alone by ilp^' ~ —ivp J^j^drTZ^ {7^)drTZ{r), 
while the triplet-superconducting term acquires the form 
H^^^ « - A ^ dr sgn{r)[TZ^ {r)TZ^-r) + h.c.]. The 
noninteracting case can therefore be written as Hq = 
J^^dr Rt(r)-HR(r), with 



n = 



-i^dr -Asgn(r) 
— Asgn(r) i^dr 



(4) 



and R(r) = [7^(r), 7^t (_r)]^. Using R(r) ^ 
(^g,'>'3Tr/'i:/y/2)J2^[ue{r),Ve{r)]'^je, where the normalized 
functions u^{r) and We(r) satisfy the eigenvalue equa- 
tion 7i[u^{r),Ve{r)]'^ = e[u^{r),v^{r)]'^ , we obtain Hq = 
^7e7e- e = there exists a localized mode at each 
edge. At r = it is of the form u^=o{r) oc g-^^l''!/^^, 
with vo{r) = mo(f). The operator corresponding to the 
edge mode, •yo ~ J druQ{r)TZ{r), satisfies the Majorana 
condition 70 = 7g. Thus the Majorana edge mode ob- 
tained by combining the right and left modes is given 
by, 



*^£oW =^^7osin(fcf^r)e-''/«, 



(5) 



for L ^ ^, where C is the normalization constant and 
^ = vf/2A the localization length. Note that in ID the 
decay is purely exponential. It is interesting to note that 
in exact analogy with the discrete lattice model, out of 
the two possible Majorana states that can be constructed 
from the fermion field, the localized Majorana that we 
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obtain at one edge corresponds to the choice ^' + 
Analogously, the Majorana localized at the other edge 
corresponds to — "S/^/i. Moreover, similar to the edge 
modes in the discrete model, those obtained in the contin- 
uum limit vanish at alternate sites for half-filling. How- 
ever, the result obtained in the continuum limit is valid 
even away from half-filling and so more general. 

Interaction effects. Next we include interactions be- 
tween the fermions given by J drdr'V{r — r')p[r)p{r') 
with V{r) being the repulsive potential and p{r) the 
fermion density. Interactions in general reduce A, and 
as a consequence ^ increases. To analyze this effect, we 
bosonize the Hamiltonian taking into consideration that 
the low-energy physics is described by a single species of 
fermions in the subband. Using the standard proce- 
dure [3, the bosonic Hamiltonian reads. 



H 



dr 
T 



VFK{dref 



K 



4A 



U 



cos[4:\fn(f) — Akpr) 



sin(2V7r6') 
(6) 



where a is the lattice constant, the dr4> field describes the 
density fluctuations and the 6 is the conjugated field. The 
quadratic part in Eq. (j6]) includes repulsive interaction 
between the fermions {K < 1), the sine term is due to the 

('2') 

triplet superconducting term Hq given in Eq. ([3]), and 
the cosine term describes umklapp scattering by V{r). 
The umklapp terms play a role only in lattice systems but 
are absent in quasi-one-dimensional quantum wires fab- 
ricated on a two-dimensional electron gas. For fermions 
on a lattice near half-filling, 4(fcF — 7r/2a)L ^ 1 and 
the oscillatory part inside the cosine term can be ne- 
glected. The interactions then lead to the renormaliza- 
tion of the coupling constants A, U, and K, which by 
standard renormalization group (RG) theory 
pressed by the RG equations 



18| 



IS ex- 



dl 2K ^ 
dS , 1 . dy 

-di^^^-R^^' m 



(2 - AK)y, 



(7) 
(8) 



where / = In [a/ag] is the flow parameter with oq the 
initial value of the lattice constant. The dimensionless 
superconducting term at the length scale a is deflned as 
5{l) = 4aA{l)/vF and y{l) = U{l)a/'KVF- The initial 
values of the rescaled parameters are given by Kq, Aq, 
do, Uq, and j/o- For K < 1/2 the umklapp term is rele- 
vant and superconductivity irrelevant, leading to a Mott 
phase, whereas for K > 1/2 the opposite is true and the 
system is superconducting. Near K = 1/2 the low-energy 
physics depends critically on the relative strength of 6a 
and 2/0- A large So compared to yo favors superconduc- 
tivity over the Mott phase and vice- versa. An interesting 
scenario corresponds to the line of flxed points Sq ~ ya 
and Kq = 1/2, where the parameters remain invariant 



under the RG flow. Following Refs. [l^, [l^, we find that 
under a change of quantization axis the theory is de- 
scribed by a quadratic Hamiltonian. Therefore, similar 
to the discrete model with t = A ~ U /4, the spectrum 
is gapless. The Majorana edge states are thus absent on 
the line of fixed points, as well as in the Mott phase. 

Away from half-filling, the umklapp term in Eq. (|6]) 
becomes strongly oscillating and can be neglected, al- 
lowing us to set y = in Eq. (|8]). The remaining 
RG equations reduce to the standard Kosterlitz-Thouless 
(KT) equations under the change of variables K 1/2^" 
and S d/V2 The flow equation of A{1) is, 

dA/dl = (1 — A'~^)A, and its solution in terms of K{1) 
is given by 



A V8[^(0 - A-q] - 41n[i^(0/A-o] +^ 

where for small deviations of K from its initial value ATg, 
/ is given by. 



I 



— = cot 

Ja 



X\/a 



(10) 



where x ~ {K~Kq)/ Kq, ko — 2A'o— 1, and a = 6^/2^ 
Rather than linearizing the KT flow eqs. around the 
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FIG. 1: RG flow of A/Ao as a function of K for Ao = 

0. 051;F/fflo and the three initial values Ko = 0.5, Ka = 0.6, 
and Ao = 0.8. Tlie solid lines are obtained from the numeri- 
cal integration of the KT eqs. The dashed lines are obtained 
from Eqs. ((9]) and (|10p [the dashed hne with the steepest de- 
cay for Ao = 0.5 is obtained from Eq. @ and I ~ {2Ko/Sq)x 

1 . The flow reaches the non-interacting hmit at A = 1 (shown 
by the red dotted line). The vertical arrows indicate the po- 
sition where 5 = 1 is reached. 



flxed point as is often done [l8|, the solutions given by 
Eqs. (O and pUj) are obtained by integrating the KT 
equations. Figure [T] shows A/Aq as a function of K 
for Aq = 0.05vF/ao and three different values of Kq, 
Kq = 0.5, 0.6 and 0.8. For all the A'o's considered, A 
reduces from its initial value and acquires its minimum 
at AT = 1. Note that near AT = 1, A shows very little 
variation. For the strongly repulsive case, Kq = 0.5, A is 
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reduced by an order of magnitude as K reaches JiT < 1 . 
In particular, for K « 0.5 and x ^ 1, Eq. pO)) can be 
approximated as I « {2Kq/Sq)x and thus A has an ex- 
ponential drop. More generally, the exponential decay 
persists as long as a; ^ ^o/(2 max{fco, a/H}) is satisfied. 
At a; ^ (5q/(2 maxjfco, -\/H}), one has to consider the 
full form for I as given by Eq. ([T0| . 

Next we discuss in detail the RG flow of the param- 
eters and its consequence for the Majorana edge states. 
Although everywhere in the repulsive regime (K < 1) K 
has a monotonic increase and A a monotonic decrease, 
the flow can be divided into two regions based on the ini- 
tial values of 6o and Kq . The most favorable scenario for 
the existence of the Majorana edge modes corresponds to 
the initial value {Ko,6q) with Kq > 1/2 in the screened 
regime or, for Kq < 1/2 with Sq > 2y/2Ko ~ ln(2/'Sroe) 
(i.e., above the separatrix). In these regions the flow is 
towards the strong coupling regime, and although A de- 
creases monotonically it remains finite. The minimum 
is reached at the length scale a{li), where K{li) = 1, 
beyond which point A increases. The RG flow crosses 
K = 1 ii the length scale a{li) is shorter than any cut-off 
length, i.e., a{li) < min{L, Lt-, a{ls)} [where Is is de- 
fined as S{ls) ~ 1 and Lt = vp/ksT is the thermal 
length]. We note that A' = 1 is a special line where 
the interaction has scaled down to zero and the solu- 
tion is obtained exactly as in the non-interacting case 
without further resorting to the RG. The Majorana edge 
state has the same form as in Eq. (O, albeit A is now 
given by the reduced value A(/i). However, for preserv- 
ing the Majorana property, which is of particular inter- 
est for the quantum computational use of the Majorana 
edge states 0, , the two edge states must have mini- 
mal overlap, i.e., x = 2A{Ii)L/vf 3> 1. Thus the drop 
in A due to the interactions should be compensated by 
increasing the length of the wire by at least a factor of 
Ao/A(/i), where A(^i) can be evaluated from Eqs. ^ 
and (fTO| . If, however, a{l*) = min{L,LT,a{ls)} < a{li) 
then the RG will be cut-off before K ^ 1 is reached. 
In the scenario when K{1*) < 1, we note from Fig. [T] 
that A(l*) K, A(li), thus we expect that the Majo- 
rana edge state will still be described by Eq. (O with 
A = A(/*). The second regime is the unscreened regime 
with ifo < 1/2 and < 2^/27^0 - ln(2ifoe). Here the 
flow is towards the line of Luttinger-liquid fixed points, 
A = and Kq < K < 1/2. In a realistic scenario the 
flow is stopped before the fixed points are reached at a 
length scale given by, a{l*) = min{L, Lt}- If a{l*) = Lt, 
then A(l*) < ksT and thermal fluctuations overcome su- 
perconductivity. On the other hand, if the wire-length L 
is the cut-off, then the superconducting term is renor- 
malized down to A(/*) sa Ao(L/ao)^-i/^o. In either 
case the bulk spectrum remain gapless and all correla- 
tions exhibit power-law decay. Thus, the Majorana edge 
states which require the presence of gapped bulk modes 
are absent. One way to ensure a gapped phase in the 



bulk is to consider a larger value for 5q. A large Sq will 
be difficult to achieve as the proximity induced gap A5 is 
further suppressed by the small ratio, aukp/ Az- More- 
over, in contrast to Kq, controlling and scaling up the 
strength of the superconducting order parameter is non- 
trivial. A simpler alternative would be to apply gates on 
top of the wire to screen the interactions and to increase 
Kq to a larger Kq that pushes the initial point {Kq,6o) 
above the separatrix, 6q > 2^J2K'q — ln(2/\Qe) or beyond 
K'q > 1/2, so that the flow is towards the strong coupling 
regime. 

Potential candidate systems for the observation of Ma- 
jorana edge states arc the helical conductors formed at 
the boundaries of topological insulators [2^ '2l'|, InAs 
nanowircs with strong spin-orbit interaction [2. 6, 22, 23i], 
quasi- ID unconventional superconductors [9|, carbon 
nanotubes Q, and quantum wires with nuclear spin or- 
dering 2J|. The latter two systems may be particularly 
interesting because they are readily available and support 
helical modes without external magnetic flelds. 
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